Improved search for D^-D^ mixing using semileptonic decays at Belle 

U. Bitenc/2 I. Adachi/ H. Aihara,^'' K. Arinstein/ T. Aushev,!*'-" A. M. Bakich,^^ V. Balagura," I. Bedny,i 
K. Belous/'' A. Bondar,i A. Bozck,^^ M. Bracko,!^'!^ Brodzicka/ T. E. Browder,^ M.-C. Chang,^ A. Chen,^^ 
W. T. Chen,22 B. G. Cheon,^ R. Chistov," I.-S. Cho,"^ Y. Choi,^^ J. Dalseno,!^ M. Danilov," M. Dash,^^ 
A. Drutskoy,2 S. Eidclman,i S. Fratina/^ B. Golob,!^'!^ jj. Ha/" J. Haba/ T. Hara.^o K. Hayasaka,20 
H. Hayashii,2i M. Hazumi/ D. Hcffernan,^" Y. Hoshi,^^ W.-S. Hou,24 Y. B. Hsiung,24 H. J. Hyun/^ K. Inami,^^ 

A. Ishikawa,32 H. Ishino,-*i R. Itoh/ M. Iwasaki,4° D. H. Kah,^^ H. Kaji,^^ P. Kapusta,^^ N. Katayama/ 
T. Kawasaki,27 H. Kichimi/ H. J. Kim,^^ Y. J. Kim,* K. Kinoshita,^ S. Korpar,i8.i2 Y. Kozakai,20 P. Krizan,!^'!^ 
P. Krokovny/ R. Kumar,^! C. C. Kuo,^^ Y. Kuroki,30 Y.-J. Kwon,*^ J. S. Lee,^^ M. J. Lee,^* S. E. Lee,^* 

T. Lesiak,25 J. Li,^ C. Liu,^^ D. Liventsev," F. Mandl,^ A. Matyja,^^ S. McOnie,^^ T. Medvedeva," 
H. Miyake,^'' H. Miyata,^^ Y. Miyazaki,^" R. Mizuk," M. NakaoJ H. Nakazawa,22 S. Nishida/ O. Nitoh,''^ 
T. Nozaki/ S. Ogawa,^^ T. Ohshima,^" S. Okuno,!^ S. L. Olsen,^'^ H. Ozaki/ P. Pakhlov," G. Pakhlova," 
H. Palka,25 C. W. Park,^^ H. Park/^ L. S. Peak,^^ R. Pestotnik,!^ L. E. Piilonen,** H. Sahoo,^ Y. SakaiJ 
O. Schneider,!^ J. Schumann/ A. J. Schwartz,^ K. Senyo,^^ M. E. Sevior,!'^ M. Shapkin,!" H. Shibuya,^'^ 
J.-G. Shiu,24 J. B. Singh,3i ^ Sokolov,io A. Somov,^ S. Stanic,^* M. Staric,!^ T. Sumiyoshi,*^ F. TakasakiJ 
N. Tamura,27 M. Tanaka/ G. N. Taylor,!^ Y. Teramoto,29 I. Tikhomirov," K. TrabelsiJ S. UeharaJ K. Veno,'^* 
T. Uglov," Y. Unno,^ S. Uno/ P. Urquijo/^ G. Varner,^ K. E. Varvell,^^ K. Vervink,!^ S. Villa,!*' 

A. Vinokurova,! C. C. Wang,^* C. H. Wang,23 M.-Z. Wang,^* p. Wang,« Y. Watanabe,!^ R. Wedd,!^ E. Won,!* 

B. D. Yabsley,^*' H. Yamamoto,^^ Y. Yamashita,^^ C. Z. Yuan,^ C. C. Zliang,^ Z. P. Zliang,33 and A. Zupanc^^ 

(The Belle Collaboration) 

' Budker Institute of Nuclear Physics, Novosibirsk 
^ University of Cincinnati, Cincinnati, Ohio 45221 
^Department of Physics, Fu Jen Catholic University, Taipei 
^ The Graduate University for Advanced Studies, Hayama 
^Hanyang University, Seoul 
''University of Hawaii, Honolulu, Hawaii 96822 
''High Energy Accelerator Research Organization (KEK), Tsukuba 
^Institute of High Energy Physics, Chinese Academy of Sciences, Beijing 
'^Institute of High Energy Physics, Vienna 
^'^ Institute of High Energy Physics, Protvino 
^''institute for Theoretical and Experimental Physics, Moscow 
^^,1. Stefan Institute, Ljubljana 
Kanagawa University, Yokohama 
^^Korea University, Seoul 
Kyungpook National University, Taegu 
^^Ecole Polytechnique Federale de Lausanne (EPFL), Lausanne 
^''Faculty of Mathematics and Physics, University of Ljubljana, Ljubljana 
University of Maribor, Maribor 
University of Melbourne, School of Physics, Victoria 3010 
^"Nagoya University, Nagoya 
Nara Women's University, Nara 
National Central University, Chung-li 
'^"'National United University, Miao Li 
''^'^ Department of Physics, National Taiwan University, Taipei 
^^H. Niewodniczanski Institute of Nuclear Physics, Krakow 
^''Nippon Dental University, Niigata 
'''Niigata University, Niigata 
University of Nova Gorica, Nova Gorica 
Osaka City University, Osaka 

Osaka University, Osaka 
Panjab University, Chandigarh 
''^^Saga University, Saga 
University of Science and Technology of China, Hefei 
Seoul National University, Seoul 
'^'^ Sungkyunkwan University, Suwon 
University of Sydney, Sydney, New South Wales 
'^''Toho University, Funabashi 
^^Tohoku Gakuin University, Tagajo 



2 



Tohoku University, Sendai 
^"Department of Physics, University of Tokyo, Tokyo 
Tokyo Institute of Technology, Tokyo 
^'"^ Tokyo Metropolitan University, Tokyo 
^'^ Tokyo University of Agriculture and Technology, Tokyo 
Virginia Polytechnic Institute and State University, Blacksburg, Virginia 24-061 

Yonsei University, Seoul 
(Dated: April 22, 2008) 

A search for mixing in the neutral D meson system has been performed using semileptonic 
K^*)~e^i/ and — > K^*^~ fj.'^u decays. Neutral D mesons from D*^ —* D^vr^ decays are used and 
the flavor at production is tagged by the charge of the slow pion. The measurement is performed 
using 492 fb~^ of data recorded by the Belle detector. From the yield of right-sign and wrong-sign 
decays arising from non- mixed and mixed events, respectively, we measure the ratio of the time- 
integrated mixing rate to the unmixed rate to be Rm ~ (1.3 ± 2.2 ± 2.0) x 10"''. This corresponds 
to an upper limit of Rm < 6.1 x lO"** at the 90% C.L. 

PACS numbers: 14.40.Lb, 13.20.Fc, 12.15.Ff 



I. INTRODUCTION 

The phenomenon of mixing has been well established 
in the K°-K°, B°-W and B°-^s° systems. Recently, 
evidence for mixing in the D^-D'^ system has been ob- 
tained with a statistical significance of more than three 
standard deviations for the first time [J, 01 ■ In addition, 
several new measurements help constrain the relevant 
mixing parameters |3| • The parameters used to char- 
acterize D^-D^ mixing are x = Ato/F and y — AF/2F, 
where Am and AF are the differences in mass and de- 
cay width between the two neutral charmed meson mass 
eigenstates, and F is the mean decay width. The mixing 
rate within the Standard Model is expected to be small 
: the largest predicted values, which include the impact 
of long distance dynamics, are of the order \y\ ^ 10^^. 

For x,y ^ 1 and negligible CP violation, the time- 
dependent mixing probability for semileptonic Z)" decays 
has the following form 



X+rVe) oc Rm 



-rt 



(1) 



where Rm is the ratio of the time-integrated mixing prob- 
ability to the time-integrated non-mixing probability: 



Rm 



dt ViDO X-l+vi) 



(2) 



The mixing rate Rm can be measured directly by using 
semileptonic decays of D° mesons. The most stringent 
constraint from semileptonic decays, Rm < 1-0 x 10""^ at 
the 90% confidence level, comes from our previous mea- 
surement i7i]. Other measurements oi Rm using semilep- 
tonic decays are less sensitive [1, H, [Tol . whereas results 
from hadronic decays are more precise[lll [T^ . IT3| . In this 
paper we present an improved search for D^-D'^ mixing 
using semileptonic decays of charmed mesons, which su- 
persedes our previous measurement 7:]. We measure Rm 
in a 492 fb~^ data sample recorded by the Belle detec- 
tor at the KEKB asymmetric-energy e+e~ collider [3], 
at a center-of-mass (cms) energy of 10.58 GeV. The Belle 7r+, and £+, where £^ can be either a muon or an 



detector [15| is a large-solid-angle magnetic spectrometer 
that consists of a silicon vertex detector (SVD) , a 50-layer 
central drift chamber (CDC), an array of aerogel thresh- 
old Chcrenkov counters (ACC), a barrel- like arrangement 
of timc-of-flight scintillation counters (TOF), and an 
electromagnetic calorimeter (ECL) comprised of CsI(Tl) 
crystals located inside a superconducting solenoid coil 
that provides a 1.5T magnetic field. An iron flux-return 
located outside of the coil is instrumented to detect Kl 
mesons and to identify muons (KLM). Two different in- 
ner detector configurations were used. The first 140 fb~^ 
of data were taken using a 2.0 cm radius beam-pipe and 
a 3-layer silicon vertex detector (SVD-1), and the subse- 
quent 352 fb^^ were taken using a 1.5 cm radius beam- 
pipe, a 4-layer silicon detector (SVD-2) and a small-cell 
inner drift chamber p^ . 

To study signal and background distributions we use 
Monte Carlo (MC) simulated samples 17] in which the 
number of selected events is about 2.7 times larger than 
in the data sample. 



II. RECONSTRUCTION OF DECAYS 

We select D'^ mesons arising from D*+ D^T^t de- 
cays and reconstruct them as D'^ K~t^vi, where 
can be either an electron or muon [l^ . The notation 7r+ 
denotes a slow pion, i.e., the pion that originates from the 
The average momentum of this pion is only about 
0.23GeV/c, whereas the average momentum of the lep- 
ton and kaon from the signal decay are 0.96 GeV/ c and 
1.52GeV/c, respectively. The momenta given in this pa- 
per are measured in the laboratory frame, unless oth- 
erwise stated; momenta measured in the cms frame are 
denoted with an asterisk, e.g., p* . The reconstruction of 
mesons in this specific decay chain enables tagging of 
the D'^ / meson fiavor at production using the charge 
of the slow pion tt^ . 

There are three detected particles in the final state: 
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electron. The non-mixed decay results in a charge com- 
bination 7r+ K~ , which we refer to as the Right-Sign 
(RS) charge combination. The mixing process results in 
7r+ £~, which we refer to as the Wrong-Sign (WS) 
charge combination, as summarized in Table |T1 

TABLE I: The definition of the Right-Sign (RS) and Wrong- 
Sign (WS) charge combinations. 



charge 


combination 


process name 






non-mixed Right-Sign, RS 




K+, £- 


mixed Wrong-Sign, WS 



Because the neutrino is not directly reconstructed, 
the masses of the D° and D*'^ candidates are smeared. 
However, by calculating the difference between the two 
masses, the uncertainty due to the neutrino four momen- 
tum cancels to a large extent. Thus 

AM = M{tTsK£v) - M{K£v), (3) 

the reconstructed invariant mass difference between the 
D*"*" and the meson, is the most appropriate observ- 
able to extract the number of signal events. For signal 
events, the distribution of AM peaks at 0.145 GeV/c^, 
the mass difference between the Z?*^ and meson (see 
Fig.[T]). 

A. Selection criteria 

Among all the different processes occurring in e+e^ 
collisions, hadronic final states are selected with an effi- 
ciency above 99%. The selection is based on the energy 
of the charged tracks and neutral clusters, total visible 
energy in the cms system, the z component (opposite to 
the positron beam direction) of the total cms momen- 
tum, and the position of the reconstructed event vertex 

Using MC simulation, the criteria to select the signal 
decays are optimized to give the best significance for the 
extracted number of mixed (WS) events, TV-^g/cr^bkg. 
The uncertainty u^-bkg, is due to the fluctuation of the 

background in the region AM < 0.16GeV/c^; the fluc- 
tuation of the signal events at the rate of our previously 
measured upper limit 7] is negligible. Hence we maxi- 
mize 



/,rbkg,AA^<0.16 

Since the kinematic properties of mixed and non-mixed 
events are the same, in the optimization the RS signal is 
used instead of the WS signal. The optimal values for 
the selection criteria in some observables are correlated 
and hence the final criteria are obtained by iterative opti- 
mization. In the optimization, selection criteria based on 



the proper decay time are also included, as described 
in Sec. [ml 

We suppress mesons arising from T(4S') — ^ BB 
events in order to avoid the situation in which the se- 
lected sample would be composed of two subsamples with 
different kinematic properties. These BB events have dif- 
ferent kinematic properties from the decays of mesons 
produced in e^e~ — > cc (continuum events), and a differ- 
ent apparent decay length between the interaction point 
and the D'^ decay vertex, because of the finite B life- 
time. As a result, the mesons from this source have 
slightly different resolutions in kinematic variables, and 
their proper decay time cannot be measured in the same 
way as for D'^ mesons from the continuum. Since the 
fraction of mesons from B decays is smaller than that 
from the continuum production, and the background con- 
tribution from B decays is large, the sensitivity to mixed 
events is not reduced by rejecting candidates from B de- 
cays. 

The quantity used to discriminate between BB events 
(spherical) and continuum events (jet-like) is the ratio of 
the second to zeroth Fox- Wolfram moment, i?2 [l^l- To 
suppress candidates from B decays we demand i?2 > 0.2. 
A further effective rejection of BB events is described 
below. 

Tracks with an impact parameter with respect to the 
interaction point in the radial direction, dr < 1 cm, and 
in the beam direction, \dz\ < 2 cm, are considered as nf 
candidates. These criteria remove badly reconstructed 
tracks and tracks not arising from the interaction point. 
A slow pion candidate is required to have a momentum 
smaller than 600MeV/c. To reduce the background from 
electrons, we require the electron identification likelihood 
(based on the information from the CDC, ACC and ECL 
|2l| ) of a TTs candidate to be £e < 0.1, which selects 
slow pions with an efficiency of 96% and rejects 72% of 
electrons. The total efficiency of the slow pion selection 
criteria and tracking is 51%. 

Electron candidates are required to have momenta 
greater than 250MeV/c and an electron identification 
likelihood Ce > 0.95. The efficiency of the identifica- 
tion criterion is 76%; in total around 46% of all gen- 
erated signal electrons are retained. Muon candidates 
are required to have momentum greater than 650MeV/c 
and the muon identification likelihood (based on infor- 
mation from the KLM and properties of the track (2^ ') 
> 0.97; the latter criterion selects muons in the chosen 
momentum range with an efficiency of 67%. These two 
requirements are highly correlated since the identification 
efficiency of muons with momenta lower than 600MeV/c 
is very poor. The efficiency of the above selection criteria 
and tracking is 30%. 

Kaon candidates are required not to satisfy the lepton 
selection criteria. Kaons from _D" — ^ Kei' decays should 
have p > 850MeV/c and kaons from Kjiv decays 

p > 600MeV/c. The difference in this requirement is 
due to the correlation between lepton momenta and kaon 
momenta that enters through other kinematic variables. 
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and due to different background contributions in both 
decay modes. A combined likelihood for a given track to 
be a , tt^ or is obtained based on the information 
from the TOF, CDC and ACC [H. Kaon candidates 
are selected using c(K)+c(Tr) ^ ^-^^ (efficiency of 87% 
for signal kaons in the selected momentum range) and 



C(K)+C(p) 



> 0.01 (efficiency of 99% for signal kaons). 
Around 42% of all generated kaons in the electron decay 
mode, and around 48% in the muon decay mode pass the 
selection criteria. 

At this stage, about 17.0% of all generated 
TT+D°,D" K-e+Ue decays, and about 12.5% of aU 



generated D 



*+ 



"D^^D^ — > K ^JL^v^ decays, are re- 



constructed. Further criteria are applied to improve the 
sensitivity to mixed events. In the following, these crite- 
ria are described and for each of them the signal loss and 
the background rejection factors are given. 
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FIG. 1: AM distribution for MC-simuIated signal events 
in different p'^i bins: pj^^ < 2.0GeV/c (dashed line), 
2.0GeV/c < p*jfi < 3.2GeV/c (solid line) and p*jfi > 
3.2GeV/c (dotted line). The histograms are normalized to 
the same area. The resolution is improved at higher values of 
momentum. The left plot is for the electron decay mode and 
the right one for the muon decay mode. 
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FIG. 2: MC simulated distribution of pj^^ for signal (solid 
line) and background events {bb, cc and uds components of 
background are shown). The arrow shows the value of the 
pl^fi requirement. 

The most effective requirement is the one on the sum of 
the kaon and lepton momenta, calculated in the cms sys- 
tem, p*j^g, see Fig. [2] Its optimized value is between 1.7 



and 1.9GeV/c. However, simulated signal events show 
a clear improvement in the AM resolution at higher 
P*Ke values, see Fig. [TJ The full width at half maxi- 
mum (FWHM) of the simulated AM distribution reduces 
from 8.6MeV/c for < 2.0GeV/c to 6.9MeV/c for 
P*Ke ^ 3.2GeV/c in the electron decay mode, and from 
7.7MeV/c to 5.2 MeV/c for the same intervals in the 
muon decay mode. Hence is required to be at least 
2.0GeV/c, a value that also eliminates a large fraction 
of D'-^ meson decays arising from T(4S') BB events. 
This requirement results in a signal loss of 28% in the 
electron decay mode and 23% in the muon decay mode, 
while rejecting 76% of the total background in the elec- 
tron decay mode and 67% of the total background in the 
muon decay mode. 

We apply a selection on the invariant mass of the kaon- 
lepton system. For the electron decay mode the optimal 
range is 0.9 GeV/c^ < M{Ke) < 1.75 GeV/c^ and for the 
muon decay mode l.OGeV/c^ < M{K^) < 1.75GeV/c'^, 
see Fig [3] In the electron decay mode this requirement 
rejects 25% of the total background at a cost of losing 
5.5% of signal events. In the muon decay mode the signal 
loss is higher, 12%, but so is the background rejection, 
44%. 
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FIG. 3: The distribution of M{K£) for signal (solid line) 
and for background events (dashed), normalized to the same 
number of entries. The arrows show the values of the selection 
criteria. The left plot is for the electron decay mode and the 
right one for the muon decay mode. 

D*^ decays to two mesons in the final state are an im- 
portant source of background. There are four such de- 
cays (see Table IH)) . Their branching fractions are at least 
one order-of-magnitude larger than the effective branch- 
ing fraction for K~^£~i'e. Although particle iden- 
tification reduces their presence in the final sample, it is 
still important to suppress this background, because its 
AM distribution has a shape similar to that of the sig- 
nal. Background exhibiting a peak around 0.145 GeV/c^ 
is called peaking background and reduces the sensitivity 
to mixed events much more than the non-peaking back- 
ground. 

In the electron decay mode, the requirement on M{K£) 
rejects 90% of the background arising from the dou- 
bly Cabibbo suppressed decay D'^ — *■ J^+vr^, through 
misidentification of the pion as an electron; in the muon 
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TABLE II: Two-body decays, representing a source of a peak- 
ing background. The symbol represents a misidentifica- 
tion. 

decay mode Br flO'^l [23] contribution to WS 



D° - 



► K-TV+ 



38.0 ±0.7 R-^t 
0.143 ±0.004 TT" 
3.84 ±0.10 K' 
1.36 ±0.03 TT-^r 



K+ 




0.14 0.15 0.16 

AM [GeV/c- 



0.17 0.18 



0.1 0.2 0.3 0.4 0.5 0.6 



decay mode the suppression rate for this background is 
98%. In both decay modes it completely eliminates the 
background due to misidentification of both pions from 
Z?" 7r-7r+. 

The Cabibbo-favoured decays 13° — > K^tt^ contribute 
to the WS background if the kaon is misidentified as 
a lepton and the pion is misidentified as a kaon. To 
suppress this type of background, the invariant mass of 
the kaon-lepton system, M^/^ (X^), is calculated with the 
pion mass assigned to the kaon candidate and the kaon 
mass assigned to the lepton candidate. If \M.,^K{Ke) — 
mjjo] < lOMeV/c^ in the electron decay mode, and 
\MTrK{Kfi) — m£)o| < 15MeV/c^ in the muon decay 
mode, the K — H candidate is rejected. Here mjjo is 
the mass of the D° meson, 1.8645 GeV/c^ [11]. In the 
electron decay mode, this requirement rejects (64 ± 2)% 
of the WS background from A'~7r+ decays, and 

(85 ± 3)% in the muon decay mode. 

To suppress the contribution from D'^ K^K^ de- 
cays {K~ being misidentified as a lepton), the invari- 
ant mass of the kaon-lepton system, Mkk{K£), is cal- 
culated with the kaon mass assigned to both candidates. 
If \MKK{Ke) - moo \ < 10 MeV/c^ in the electron decay 
mode, and \Mkk{Kij.) — mjjo \ < 15 MeV/c^ in the muon 
decay mode, the K — £ candidate is rejected. In the elec- 
tron decay mode, this requirement rejects (70 ± 5)% of 
the WS background from 13° K~ decays, and in 
the muon decay mode (89 ± 1)%. 

The requirements on M.^K{Kt) and Mkk{K£) result 
in a signal loss of 3% in the electron decay mode and 
2% in the muon decay mode. The rejection of the total 
background in both decay modes is similar to the signal 
loss. 



B. Rejection of 7 ^ e e 

An important source of background is due to electrons 
from photon conversions: either the electron candidate, 
the slow pion candidate, or both, may be due to 7 — > 
e^e'' tracks. 

In the electron WS sample, both the slow pion and the 
signal electron candidates can come from 7 e+e^, and 
such events tend to have low AM values (Fig. 01 left). 
To suppress this background we calculate Mee(7rsesignai), 
the invariant mass of the tt^ — Csignai system with the elec- 



FIG. 4: Left: the AAf distribution of 7 — » e+e" events, 
where one of the electrons is selected as the slow pion candi- 
date and the other as the electron candidate, before applying 
the selection on Mee(7rsesignai). Right: the Mee(7rsesignai) dis- 
tribution for these events (solid line) , for signal events (dotted 
line) and for the total WS background (dashed line) . The ar- 
row shows the value of the selection criterion. 



tron mass assigned to both tracks. We reject candidates 
with Mee(7rsesignai) < 0.14GeV/c^. The requirement re- 
jects more than 99% of this background, see Fig.lH right. 
To retain equal reconstruction efficiencies for the mixed 
and non-mixed events, this requirement is implemented 
in both RS and WS samples. 

Background from t^/tt^ 77 — > e+e^ e+e^, where 
one electron from the first photon and one electron from 
the other photon are taken as the electron and slow pion 
candidate, exhibits similar behavior to the case where 
the electrons are both from the same photon. This back- 
ground is also successfully eliminated by the above selec- 
tion. 

Assuming that the signal electron candidate comes 
from 7 — !■ e+e^, in the electron decay mode we per- 
form a search for the other electron 62 among all the 
other tracks in the event with the opposite charge to the 
signal electron candidate. If M(esignaie2), the mass of 
the Csignai — 62 systcm, is below 80MeV/c^, the electron 
candidate is rejected. 

Assuming that the slow pion candidate is a misidenti- 
fied electron from 7 — > e+e^, we perform a search for 
the other electron (62) among all the other tracks in 
the event. The other electron should have the opposite 
charge to the slow pion candidate, and an electron like- 
lihood Ce > 0.8 to reduce rejection of true signal slow 
pions. If Mee(7rse2), the mass of the tTs — 62 system 
with the electron mass assigned to both tracks, is be- 
low 80MeV/c^, the slow pion candidate is rejected. This 
photon conversion rejection is performed for slow pion 
candidates in both the electron and muon decay modes 
and results in around 0.4% signal loss and rejects around 
2.4% of the total background. 

In total, the rejection of photon conversion in the elec- 
tron decay mode results in a 14% signal loss and 32% 
rejection of the total WS background. 
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Neutrino reconstruction 



Four-momentum conservation in e^e collision implies 



(5) 



for the signal decay, where Pcms stands for the cms 
four-momentum of the e^e~ system and Prost indicates 
the four-momentum of all detected particles except the 
charged kaon and the lepton candidates [13] ■ Eq. ([5]) is 
true if all the particles produced in the e+e" collision are 
detected. As the Belle detector covers nearly the entire 
solid angle around the interaction point, neutrino recon- 
struction can be successfully performed. 

The variable Prcst is calculated using all the remaining 
charged tracks (except the kaon and lepton candidates) 
with dr < 2 cm and \dz\ < 5 cm, and photons with an en- 
ergy above 70 MeV. Mass is assigned to a track according 
to the following criteria: 

• A track is assigned the electron mass if its electron 
likelihood is Ce > 0.9. 

• A track is assigned the muon mass if Ce < 0.9 and 
its muon likelihood is Cf^ > 0.9. 

• A track is assigned the kaon mass ii Ce < 0.9, Cf^ < 



0.9 and 



C{K) 



C{K)+C{w) 



> 0.5. 



A track is assigned the proton mass if Ce < 0.9, 



Cf, < 0.9, 



< 0.5 and 



> 0.5. 



• In all other cases the track is assigned the charged 
pion mass (2^ . 

A first approximation for the neutrino four-momentum 
is obtained using Eq. ([5]) and the resulting AM dis- 
tribution for signal events is shown in Fig. [5] (left) with 
the dashed line. It peaks at around 0.148 GeV/c^, a value 
close to the D*+ mass difference, 0.145 GeV/c^, and 
has a FWHM of SSMeV/c^. 

Two kinematic constraints are used to improve the res- 
olution on the neutrino momentum. To simplify the ex- 
pressions, we performed the calculation in the cms sys- 
tem, since pl^^ = 0. First, the squared invariant mass 
of the selected particles is calculated using M^{K£i') = 
{P* + P*i^i,f/c^. The distribution of M'^{K£i^) is shown 
in Fig. [6l left. For signal events, the invariant mass 
should equal rrijjo. To reject poorly reconstructed events, 
exhibiting a large FWHM of the final AM distribu- 
tion, only candidates with -25Gey^/c^ < M^{Kiv) < 



64GeV^/c^ are retained. For the selected events, P*g^^. is 
rescaled by a factor ^ requiring 



M^iKiiy) = (P* 



(6) 



The neutrino four-momentum is then recalculated as 
-P^* = ^cms - Pxe - iP*csu and a corrected M{-KsK£v) is 
obtained, where M{K£v) has been forced to equal to^io. 
With this correction, the AM distribution has a FWHM 
of 11 MeV/c^ in the electron decay mode and lOMeV/c^ 
in the muon decay mode; the improvement is shown in 
Fig. O (left). The distribution of the scale factor ^ for 
events in the finally selected sample is shown in the left 
plot of Fig. [71 It peaks at around 1.04 for the electron 
decay mode and 1.06 for the muon decay mode. The 
average ^ in both decay modes is around 1.3. 
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FIG. 5: Distribution of AM for signal (left) and background 
events (right): with the first approximation for the neutrino 
four-momentum (dashed), after applying the constraint on 
the mass (dotted) and with the final neutrino momentum, 
obtained as described in the text (solid line). Selection criteria 
on M'^{K(.v) and M^ have been omitted. The plot is for the 
electron decay mode; the distributions in the muon decay 
mode are similar. 



As a second kinematic constraint, the square of 
the missing mass, M^, is used. The distribution of 
M^ is shown in Fig. [HI right. For events satisfying 
-5GeVVc4 <Ml < 0.5GeVVc^ the angle a between 
the direction of p^^^^^ and the direction of is corrected 
in order to yield 



(P*)2 = (p* 



P, 



Kl 



0; 



(7) 



expressed in terms of energies and magnitudes of three- 
momenta this yields 



Mf,c 



2„4 



E 



Kt 



*2 2 



2 c2 



2pJf£^P*ostC cos a = 0. 



(8) 



The angle a is corrected by rotating p^^gt in the plane determined by the vectors p^^st and The distribu- 
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FIG. 6: Left: distribution of M^{Klu), calculated usins the 
first approximation neutrino momentum. Right: distribution 
of , calculated with the neutrino momentum obtained af- 
ter applying the constraint on M^{K£u). The solid line is 
for signal and the dotted one for the background. The distri- 
butions are shown for the electron decay mode; in the muon 
decay mode they are similar with a slightly smaller root mean 
square values. 

tion of the correction angle ofNEW — Q^old for the finally 
selected signal events is shown in the right plot of Fig. [T] 
It has a peak at around 1° and an average value of 8°. 
The final neutrino four-momentum is calculated with the 
rescaled and rotated P*^^^, using Eq. ([5|). 

The requirements on {Ktv) and result in a sig- 
nal loss of 4.5% in the electron decay mode and 4.1% in 
the muon decay mode while rejecting 9.7% in of back- 
ground in the electron decay mode and 8.9% of the back- 
ground in the muon decay mode. 

The AM distribution obtained using the neutrino 
four-momentum after the use of kinematic constraints 
is shown in Fig. [5] (left) as the solid line. The resolution 
is significantly improved, with the FWHM being about 
6.6MeV/c^ for the electron decay mode and 6.2MeV/c^ 
for the muon decay mode. Using the MC-simulated back- 
ground events, it has been verified that such a neutrino 
reconstruction does not induce any peaking in the back- 
ground AM distribution. From the right plot in Fig. [H] it 
can be seen that the number of background events in the 
signal region {AM < 0.16GeV/c'^) after applying the 
constraints only slightly exceeds the number of events 
without the constraints (by 8.1% in the electron decay 
mode and 16.6% in the muon decay mode). 

For D*+ candidates we require AM < 0.18GeV/c2, 
which retains 97.4% signal events in the electron decay 
mode and 97.8% in the muon decay mode. 



D. D combined with both tt^ and 

After all the above requirements are applied, a small 
fraction of events contain a D'^ candidate that has been 
combined with slow pions of opposite charges to form 
both RS and WS D* candidates. Events in which such 
candidates are found are rejected. This veto results in a 
signal loss of about 2.6% in both electron and muon decay 
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FIG. 7: Left: distribution of the scale factor ^ obtained by the 
constraint on M^{Ktv). Right: distribution of the rotation 
angle for p^cstj obtained by the constraint on M^. Both are 
for signal events; the solid line shows the distributions for 
electron decays and the dotted one for muon decays. 

mode. For the former it rejects 18% of background events 
and for the latter about 13%. 

At this stage the efficiency for reconstructing the signal 
with AM < 0.18GeV/c2 is found to be (8.0 ± 0.3)% for 
the electron decay mode and (7.2 ± 0.3)% for the muon 
decay mode. The quoted errors include a small variation 
of efficiency depending on the detector conditions. 

A summary of the applied selection criteria is presented 
in Table IIIII Since the non-mixed and mixed processes 
have the same kinematics, the WS and RS efficiencies for 
all the described selection criteria are the same. 



III. PROPER DECAY TIME 

As the proper decay time distribution of WS back- 
ground events tends to have lower values than that of 
WS signal events, the proper decay time of a meson 
can be used to select possible mixed events with a higher 
purity. Since the information on the proper decay time is 
used only to increase the sensitivity to WS events, model- 
ing of the proper decay time distribution is not detailed. 
We do not account for the fact that the associated sig- 
nal decays (defined in Sec. IIV[) have a slightly different 
proper decay time resolution function than the signal de- 
cays, or for the fact that the resolution depends slightly 
on the true value of the proper decay time. The effects of 
these two assumptions were studied carefully; the differ- 
ences between data and the modeling functions are taken 
into account in the systematic uncertainties and lead to 
a negligible change of the final result. 

The dimensionless proper decay time (proper decay 
time in units of r^jo = (410.1 ± 1.5)ps [1^) is calculated 
from the flight distance / and its momentum p^o : 



TDOPDO 
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TABLE III: Summary of all the applied selection criteria. 





electron decay mode 


muon decay mode 


event 


i?2 > 0.2 


TVs 


impact parameter: |5r < 1 cm, \&z\ < 2 cm 
Pn, < 0.6 GeV/c 
electron likelihood < 0.1 

Mee(7r,e2) > SOMeV/c^ 


£ 


Pe > 0.25 GeV/c 
electron likelihood > 0.95 

M(esignaie2) > 80MeV/c2 


> 0.65 GeV/c 
muon likelihood > 0.97 


K 


PK > 0.85GeV/c 

LiK)/ [£{K) 
CiK)/ [LiK) 


PK > 0.6GeV/c 
f £(7r)] > 0.51 
+ C{p)] > 0.01 


K — t 


Pki > 2 

0.9GeV/c2 < M{Ke) < l.TSGeV/c^ 
\M^K{Ke)-mi:,o\ < lOMeV/c^ 
\MKK{Ke) - moo\ < lOMeV/c^ 


GeV/c 

l.OGeV/c^ < M{K^j.) < 1.75GeV/c^ 
\M^K{Kfj) - mao \ < 15MeV/c2 
\MKK{K^l) -moo \ < 15MeV/c2 


TTs — e 


Mee(7rse,ignai) > 0.14GeV/c'' 


/ 


V recon. 


-25GeVVc'* < M{Kluf < 64GeV7c'' 
-5GeVVc'' < < 0.5 GeVVc* 




AM < 0.18GeV/c^ 
reject _D° combined with both tt^ and tt^ 
1.6 < < 9.0 



The observed txy distribution is smeared due to the ex- 
perimental resolution. As the data recorded with the 
SVD-2 configuration has a shghtly narrower resolution 
function than the data taken with the SVD-1 configura- 
tion, we perform the measurements of Rm separately for 
both subsamples. Thus we have four subsamples: the 
electron subsamples are denoted by e-1 and e-2 for the 
SVD-1 and SVD-2 configurations, respectively. Similarly 
the muon subsamples will be denoted as /i-1 and fi-2. 



A. Distribution of signal events 

We obtain the resolution function for signal events 
from the data, using the RS decays. To be able to do so, 
we first determine the shape of the txy distribution for 
RS background events, which is shown in Fig.[51 right, as 
the dashed hne. This distribution is also obtained from 
the data, as described in the following. 

The normalized function that describes the t^y distri- 
bution of RS background events is 



•^5fg = (1 - • { [fMt; r") + (1 - fjS{t)] ® [f.Lit b') + (1 - fi)L,{t; h\, h^)] ) + ft ■ L{t- bt). (12) 



The D° momentum is calculated by summing the mo- 
menta of the daughter particles. The D'^ flight distance 
is the distance between the D'^ production vertex, rprod, 
and its decay vertex, rdcc- The decay vertex is obtained 
by fitting the kaon and lepton tracks to a common ver- 
tex. The production vertex is obtained by extrapolating 
the D'^ momentum vector to the e~^e~ interaction region. 
The position and width of this region are determined over 
a large number of e~^e~ interactions for which the KEKB 
beam conditions do not change significantly. 

According to MC simulation, the resolution on the 
proper decay time is improved if the D° flight distance 
is calculated as the projection of the rdoc — ^prod vector 
on the normalized D'^ momentum vector. Since the in- 
teraction region is much narrower in the radial direction, 
we use only the radial components {x and y) to measure 
the proper decay time. The radial flight distance Ixy is 
calculated as 

_ 'dec 'prod' 'dec 'prodi \P D° ^ P I , „x 

Ixy — ; ■ ll"J 



The proper decay time is then evaluated as 

rriDalxy 



txy = , "^^"^-^ (11) 



Too 



J 
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It is composed as a sum of an exponential function 
E{t; T^) with the decay time r** and a delta function (their 
fractions are determined by the parameter /g), convolved 
with a detector resolution. The latter is phenomenologi- 
cally described by a sum of the Lorentz function L{t; b^) 
and an asymmetric Lorentz function La{t\h\,b^^) (both 
are explicitly given in the Appendix, Eq. ((28|) and (|29)) ): 
h\ and h^^ are their width parameters and f\ deter- 
mines their fractions in the resolution function. A wide 
Lorentz function is added to describe the decay times 
measured from badly reconstructed tracks (outliers); its 
width parameter is 6^ and its fraction in the total sam- 
ple The convolutions are performed numerically by 
substituting the integral with a sum. It has been veri- 
fied that the numerical accuracy is satisfactory, i.e. not 
affecting the result. 

To determine the eight free parameters of the function 
from the data, we divide the txy range into 15 intervals 
as shown in Fig. [51 left. In each of these 15 intervals we 
extract the number of RS background events, A^RSbkgi 
by performing fits to the AM distribution of RS events. 
The fit to the AAf distribution is described in detail in 
Sec. IIVBI The errors include the systematic error due 
to the finite number of MC simulated events and the 
uncertainty of the correlated background fraction in the 
total RS background. 

If we divide A^^g in one of the 15 txy intervals by 
the total number of background events A^Rs'bkg i^^^ 
over the 15 intervals), the obtained fraction is expected 
to agree with the integral of T^^^ (Eq. [T^ over that t^y 
interval. 



Hence we calculate A', 



RS.bkg Ws,bkg in aU 15 t^y inter- 



vals and determine the eight free parameters of T^^^ by a 
fit to these fractions. The value of the fitting function 
in each txy interval is calculated as the integral of ^i^g 
over that txy interval. The fractions and the result of the 
fit for the e-2 subsample are shown in Fig. [51 left plot; 
the obtained background distribution is shown with the 
dashed line in the right plot. The reduced values of 
the fits are reasonable, ranging from 0.4 to 4 for 7 degrees 
of freedom. The uncertainties on the fitted parameters 
of the RS background proper decay time distribution are 
taken into account when calculating the uncertainty of 
the result, Rm- 



1. RS signal distribution 

The proper decay time distribution for RS events is 
described by 



RS 



^A/'tot(/s-^5? + (l 



(13) 



where A/tot is the total number of RS events and is the 
signal fraction, obtained from a fit to AM in the entire 
proper decay time region. The fs values are (69.9±0.2)%, 
(70.8 ± 0.1)%, (62.7 ± 0.2)% and (62.6 ± 0.1)% for e-1, 
e-2, fi-1 and /i-2 subsamples, respectively. The function 



J^lfg describes the shape of the RS signal events, which 
is an exponential convolved with the resolution function. 



T!^^ = E{t;r')^nsig. 



(14) 



The resolution function is parameterized as T^sig = 
flLait; bt, bl) + f^G{t; a) + [l - fl - fi)Lit; b^). (15) 
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FIG. 8: Left: the fraction of RS background events in 15 
proper decay time intervals (error bars) and the result of the 
fit, described in the text (dotted lines). Note that the t^y in- 
tervals are not equidistant, which causes the apparently non- 
smooth shape of the function. The first and the last interval 
include events with tj^y < —10 and txy > 10, respectively. 
Right: the txy distribution of RS events (error bars) with the 
result of the fit (solid line) for the e-2 subsample. The dashed 
line shows the RS background distribution, J-^^, obtained 
from the data as explained in the text. 

Here G{t; a) is the Gaussian function. The six free 
parameters of the resolution function for the RS signal 
events are /f , /|, bf, b^, a and 6^. They are obtained by 
a fit of J-^^ to the proper decay time distribution of 
RS events. In this fit, the parameters of are fixed to 
the values previously obtained and r*' , the dimensionless 
D'^ decay time, is fixed to 1.0. An example of a fit re- 
sult is shown in the right plot of Fig. [51 The reduced 
values range between 1.2 and 2.4 for 194 degrees of free- 
dom and exhibit a slight disagreement between the fitting 
model and the data. The disagreement is accounted for 
in the systematic error evaluation, resulting in a negligi- 
ble change of the final result. 



2. WS signal distribution 

Since decays of the mixed and of the unmixed 
D'^ mesons have the same kinematic properties, the 
proper decay time resolution function for both is assumed 
to be the same. Hence from the RS signal resolution func- 
tion, T^sig, the proper decay time distribution for WS 



signal events is calculated: 



•WS 
sig 



At^ 



' "^sigj 



(16) 
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FIG. 9: Left: the proper decay time distribution of the WS 
data events (error bars) and the distribution for the WS sig- 
nal events (sohd hue), obtained from the data according to 
Eq. (|16p . The distributions are for the e-2 subsample. Right: 
the figure-of-merit dependence on the lower and upper limit 
of the selected t^y interval. When plotting the dependence on 
the lower limit, the upper limit was set at its optimal value 
and vice versa. The arrows show the value of the selection 
criterion, Eq. (|17p . Distributions for other three subsamples 
are similar. 



where A is the corresponding normalization constant and 
T is fixed to 1. 

To select the t^y interval with the highest sensitivity to 
mixed events, the ratio given in Eq. (0]) is maximized. In 
this optimization we use the calculated distribution for 
WS signal events, J'J^^ (Eq. dH])); for the background 
we use the distribution of the WS events from the data. 
Even if the latter contains some mixed events, the effect 
on the result of the optimization is negligible, since D'^ 
mixing is small. The distributions of background events, 
WS signal events and the figure-of-merit can be seen in 
Fig. The optimal proper decay time intervals for 

the four subsamples range between 1.6-1.7 and 8.9-9.5. 
In order to keep the measurement method uniform, we 
select a common interval for all four subsamples, 



1.6 < t^y < 9.0. 



(17) 



In this interval, about 70% of WS signal events are se- 
lected, while rejecting about 80% of background events 
(the values are similar in each of the four subsamples) . 



B. Extraction of Rm and further improvement 

Because the proper decay time distribution of the RS 
signal events (e"*/"^" ) is different from that of the WS sig- 
nal events (i^e~*/'^^), after the application of a selection 
based on proper decay time, the ratio Rm is obtained as 



'I 

RS ^WS 



(18) 



where A'ws,RS are the numbers of extracted signal events 
without the t^y selection, and j^g are the numbers of 
extracted signal events in the selected t^y interval. The 



superscript i labels different t^y intervals. The efficien- 
cies e|^g are obtained by integrating the proper decay 
time distribution of the RS signal events, ^^g, over the 
selected t^y interval. Similarly, the efficiencies e^g are 
obtained from the calculated proper decay time distribu- 
tion of the WS signal events, The ratios Ers/^ws 
are listed in Table IVIII The errors on the efficiencies 
quoted there include the uncertainty on the fraction of 
the RS correlated background (as defined in Sec. IIV A|) . 
the statistical and systematic uncertainty on the signal 
fraction in the RS sample, the statistical and systematic 
uncertainties in the parameters of the RS background 
and signal t^-y distributions and the uncertainties in the 
world averages of tdo and rriDo [l^ . The resulting errors 
on the measured parameter Rm are included in the sys- 
tematic uncertainty and are negligible. Imperfections in 
modelling the decay time distributions are included as a 
separate source of a systematic uncertainty as described 
below. 

To further exploit the proper decay time information, 
we divide the chosen t^y range (Eq. (jl7p ) into six inter- 
vals, with boundaries at 1.6, 2.0, 2.5, 3.1, 4.0, 5.6, and 
9.0. The binning is chosen so as to have approximately 
the same number of events in each interval. The mix- 
ing rate is measured in each of the six intervals and the 
measurements are expected to be consistent. Due to the 
additional proper decay time information the sensitivity 
of the final result is expected to be improved in compar- 
ison with the sensitivity of a single measurement in the 
total 1.6 < txy < 9.0 range. 



IV. SIGNAL YIELD EXTRACTION 

According to the MC simulation, the selected sample 
of RS events includes many candidates from semileptonic 
decays other than K~£'^i/, combined with the cor- 

rectly reconstructed slow pion. The most important of 
these decays are: 

• D° ^ K*-£+i^e, followed by K*-^ R-n^, 

• p^l^Vf^, followed by 7r~7r°, 

• D" ^ K*-l+ve, followed by K*- Wii- . 

The final state lepton in these decays is of the same 
charge as in the decay 13° K~i'^i/, so its charge can 
be used to tag mixing in the same way. In the last three 
decays with a 7r~ instead of the in the final state, the 
pion is misidentified as a kaon. The selected sample also 
includes candidates where the D'^ semileptonic decay is 
correctly reconstructed, but the slow pion decays in flight 
to a muon, nf — > fi^v^, and then the muon is misiden- 
tified as the slow pion. The muon has the same charge 
as the slow pion and can be used to tag the flavor at 
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production in the same way as the slow pion. Hence all 
these processes are treated as part of the signal and will 
be referred to as associated signal decays. 

According to MC simulation, the associated signal de- 
cays have a similar AM distribution to the Z?" K~t^v 
decays. Due to unreconstructed or misidcntified parti- 
cles in these final states, the FWHM of the distribution 
is larger, 12.3 MeV/c^ in the electron decay mode and 
14.3 MeV/c^ in the muon decay mode (see Fig. [TO)) . For 
the same reasons, the proper decay time distribution is 
slightly different from that for the signal. The fraction 
of the associated signal decays in the sample of all recon- 
structed signal decays can be found in Table [IVl There 
is also a small fraction of signal events (around 1%) from 
BB events. Due to the lower average momentum of D 
mesons from B decays, the AM distribution for signal 
events from B decays is slightly wider than that for signal 
events from cc events; it is similar to the AM distribution 
of the associated signal decays. 
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FIG. 10: The simulated AM distribution for the signal decays 
(solid line) and for the associated signal decays (dotted line), 
for both the (left) electron and (right) muon decay modes. 



TABLE IV: The associated signal fraction and the fraction of 
signal from B decays in the total signal, in [%], as obtained 
from MC simulation. The quoted uncertainties are statistical 
only. The fractions are shown for the entire proper decay time 
interval and for l.i 





assoc. sig. [%] 


sig. from B decays [%] 




aU t^y 1.6 - 9.0 


aU t^y 1.6 - 9.0 


e 


16.58±0.05 17.7±0.1 


1.20±0.01 1.18±0.03 




11.54±0.04 12.3±0.1 


1.07±0.01 1.04±0.03 



A. Background 

The background is divided into two categories: the 
correlated background and the uncorrelated background. 
The AAf shapes of both background components can be 
seen in Fig. [TlJ 



Correlated background is the background where the 
lepton candidate or the kaon candidate, or both, origi- 
nate from the same decay chain as the slow pion candi- 
date. The angular correlation between the slow pion and 
candidates leads to a concentration of events at low 
values of AM. 

The remaining, uncorrelated background has a AM 
distribution that rises steadily from threshold, as the 
available phase space increases. This component is dom- 
inant, especially in the WS sample (see Table W\ . The 
fraction in the muon decay mode is larger than in the 
electron decay mode, because the probability for a kaon 
or a pion to be misidcntified as a muon is larger than 
the probability to be misidentified as an electron. The 
fraction is larger in the RS sample than in the WS sam- 
ple due to the larger branching fractions of the Cabibbo 
favored decays. Selecting the proper decay time interval 
1.6 < txy < 9.0, decreases the fraction of the correlated 
background in the total background. According to MC 
simulation, the correlated background has three compo- 
nents: background from Z3*+ — > tt+D" decays, which 
has the largest fraction and is the most strongly peaked 
of the three, background from — > 7!)°, 7 — » e+e" 
decays in which one of the two electrons from 7 conver- 
sion is taken as a slow pion candidate, and background 
from Kg — » TT+vr^ decays, where one of the pions is taken 
as the slow pion candidate and the other is assigned as a 
kaon or lepton candidate. 



TABLE V: The fraction of the correlated background in the 
total background in [%], as obtained from MC simulation, for 
the electron and muon decay mode, in both RS and WS sam- 
ples. Fractions for the entire proper decay time interval and 
for the selected interval are shown. The quoted uncertainties 
are statistical only. 



Right Sign sample 



aU tx 



1.6-9.0 



34.49±0.07 25.33±0.15 
40.07±0.06 39.62±0.14 



Wrong Sign sample 



all tx 



1.6 - 9.0 



7.22±0.05 5.64±0.10 
14.87±0.06 14.76±0.14 



The AM distribution for the total RS background can 
be seen in Fig. [T^] as the dashed line. To fit the AM dis- 
tribution of the data as explained in Sec. IIVBI the 
AM distribution of the RS background events is obtained 
from MC simulation. 

The AM distribution for the total WS background is 
shown in Fig. [TT] as the solid histogram. One can see 
that the difference in shape between the uncorrelated 
background and the total background is smaller in the 
electron decay mode than in the muon decay mode. This 
is both due to the larger correlated background compo- 
nent in the muon decay mode (see Table |V|) and due to 
its shape, shown in the same figure. 

To check the MC simulation and to avoid systematic 
errors arising from any discrepancy with the data, the un- 
correlated background in the WS sample was described 
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FIG. 11: The simulated AM distribution for WS background 
events. The solid line shows the total background, the dotted 
line shows the uncorrelated background and the dashed one 
shows the correlated background, multiplied by a factor of 15 
and 5 in the electron and muon samples, respectively. The 
plots are for the selected proper decay time interval. The left 
plot is for the electron decay mode and the right one for the 
muon decay mode. 



using the data. This background is modeled by combin- 
ing slow pion candidates and candidates from differ- 
ent events. Technically this is done by embedding slow 
pion candidates into other events according to the fol- 
lowing procedure: 

• All slow pion candidates from an event in which a 
D° candidate was found, are taken to be embedded 
into other events. 

• We denote by Af the number of slow pion candi- 
dates which form a WS combination in their orig- 
inal event, and their charge by Q. Slow pions are 
embedded only into events with the same value of 
A/" and only the embedded slow pions of charge Q 
are used to form the WS combinations. 

• For each value of Af, slow pion candidates from 
Ma different events are stored to be embedded into 
other events. 

• Each D'^ candidate is combined with slow pions 
from several other events; slow pions from a max- 
imum of A/a events are used. Once a combination 
satisfying all the D*~^ requirements is obtained (in- 
cluding AAf < 0.18 GeV/c^), further combinations 
are not formed. 

• Once an embedded slow pion is used to form a com- 
bination satisfying all the Z3*+ requirements, none 
of the remaining slow pions from the same event is 
embedded into other events. 

With these requirements, the AAf distribution of the 
sample of embedded slow pions slightly depends on A/a, 
the maximum number of events from which the embed- 
ded slow pions are taken and tested with a single D'^ can- 
didate. The dependence is due to the AM dependence 



on the slow pion momentum. Slow pions with higher mo- 
menta tend to form D*+ candidates with slightly higher 
AM values, hence their probability to form a D*^ can- 
didate with AAf < 0.18 GeV/c^ is smaller. Increasing 
AfA enables these slow pions to be tested with a larger 
number of D'^ candidates and enhances the probability to 
form a combination with AAf < 0.18 GeV/c^. Thus in- 
creasing A/a slightly enhances the contribution at higher 
AAf values. 

The most appropriate A/a value is determined from the 
data by observing samples of txy < 0.0. In this proper de- 
cay time region, the expected fraction of the mixed signal 
events in the WS background is much smaller (around 18 
times in the electron and around 27 times in the muon 
decay mode) than in the region 1.6 < t^y < 9.0. Hence 
it is safe to assume that the WS data sample of t^y < 0.0 
contains no mixed signal events. For txy < 0.0 we com- 
bine the SVD-1 and SVD-2 subsamples and compare the 
AAf distribution of the WS data with the AAf distri- 
bution of the subsample used to describe the WS back- 
ground, i.e. embedded slow pions with the addition of 
the MC correlated background events (as explained at 
the end of this section) . The AAf distributions are com- 
pared by observing the value of 



f 155 = A^A 



M<0.155 



/Nam 



<0.18, 



(19) 



the ratio of the number of events with AAf < 
0.155 GeV/c^ and the number of events with AAf < 
0.18 GeV/c^. This is a representative observable that 
is used to characterize the AAf distribution by a single 
number. The values for the data and for the background 
with three different values of A/a are shown in Table IVII 



TABLE VI: Comparison of the riss values for the WS data 
and the modeled WS background with txy < 0.0, using dif- 
ferent values of A/a- 



txy < 0.0 


electron mode 


muon mode 




ATa 


ri55[%] 


ATa 


ri55[%] 


data 




25.66 ±0.18 




29.00 ±0.16 


background 


40 


25.80 ± 0.05 


20 


29.11 ±0.05 




45 


25.70 ± 0.05 


25 


29.00 ± 0.05 




50 


25.61 ±0.05 


30 


28.90 ± 0.05 



From Table IVIl one can see that for the electron decay 
mode the best agreement between the data and the back- 
ground {txy < 0.0) is for TVa = 45 and in the muon decay 
mode for A/a = 25. The dependence of the final result 
on the A/a value is taken into account when evaluating 
the systematic uncertainty. 

As a cross-check, we combine SVD-1 and SVD-2 sub- 
samples and compare the AAf distribution of the MG 
uncorrelated WS background with that of the embedded 
slow pions. Their agreement is good. 

To obtain the final AAf distribution of WS background 
events, the AAf shape of the WS-correlated background 
is taken from MG simulation and added to the sample of 
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the embedded slow pions in the same fraction as found by 
MC simulation. The uncertainty on this fraction is taken 
into account when evaluating the systematic uncertainty. 



B. Fit to AM distribution 

To extract the signal yield, we perform a binned max- 
imum likelihood fit to the AM distribution, assuming a 
Poisson distribution of events in AM bins and thus max- 
imizing 



""l" r RS electron mode, SVD-2 

12 - 



12 



(a^(am,)) 



(20) 



Here Nj is the number of entries in the j-th bin and 
fi(AMj) is the expected number of events in this bin, 
given by 

m(AMj) = Mtz [fsPsiAM,) + (1 - /,)Pb(AMj)] . (21) 

Ps is the signal AM distribution obtained from MC simu- 
lation. Pb is the background AM distribution composed 
as described above. The signal fraction fs is the only 
free parameter in the fit. Afn is the number of entries 
in the fitted histogram. A'bin = 45 is the number of in- 
tervals in the AM distribution. The quoted values 

are obtained using — J2f=T ^^^^^'^^^"^^^^ , where aj 

includes the statistical uncertainties of the fitting his- 



tograms, 



C. The RS signal yield 

The fit to the AM distribution in the RS sample is 
performed as described above; examples of the fit result 
are shown in Fig. [121 In the total t^y range, the signal 
fraction fs is about 70% in the electron decay mode and 
about 63% in the muon decay mode. The fraction is 
largest for 1.6 < t^y < 2.0 (82% in the electron decay 
mode and 74% in the muon decay mode) and decreases 
at larger t^y values: for 5.6 < t^y < 9.0 it is 62% in 
the electron decay mode and 54% in the muon decay 
mode. The values of the fits in the individual t^y 
intervals are in good agreement with the expectation for 
40 degrees of freedom. In the total t^y range, the reduced 
values are larger than expected (values of 1.5-2.6 for 
40 degrees of freedom) . This is explained by a difference 
in the amount of associated signal between the data and 
the MC simulation. Repeating the fits with a fraction of 
the associated signal as the second free parameter yields 
reduced values around 1.0 also for the total t^y region. 
This effect is considered in the estimate of the systematic 
uncertainty. The numbers of RS signal events are given 
in Table ED 




RS muon mode, SVD-2 
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1 +^++ data 




1 fit 




1 backgr . 
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FIG. 12: The AM distribution of the RS events for 1.6 < 
txy < 9.0, SVD-2. The dashed line represents the background, 
the solid line is the result of the fit, described in the text, and 
the points with error bars are the data. The left plot is for 
the e-2 subsample, the right one for the /i-2 subsample. 



V. RESULT 

As the kinematic properties of the RS and WS de- 
cays are the same, we use the AM shape of the MC 
simulated RS signal events also for the WS signal de- 
cays. The AM distribution of the WS background events 
is obtained as described in Sec. IIV Al By fitting the 
AM distribution as described in Sec. IIVBI we extract 
the number of the mixed signal events in the four WS 
subsamples (e-1, e-2, /i-1, fi-2). The AM distributions 
and the x^ values of the fits for all the subsamples and 
different proper decay time intervals are shown in Figs. 
[131 [m The extracted WS signal yields are given in Table 

Ivnl 

For each of the four subsamples we determine the mix- 
ing ratio Rm by three different methods which arc dis- 
cussed below. 

1) The fit to AM in the RS and WS samples is per- 
formed without any selection based on the proper 
decay time measurement. The ratio Rm is calcu- 
lated as the ratio of the obtained number of WS 
and RS signal events, A^ws/-^RS- The results can 
be found in Table WJJ\ in rows labeled "all txy" ■ 



2) The fit to AM distributions for the RS and WS 
sample is performed for events with 1.6 < t^y < 9.0. 
The ratio Rm is calculated as N^^/N^g x ej^s/^ws- 
The results are given in Table IVIII in rows labeled 
"1.6-9.0". The resulting statistical uncertainty of 
the result is around 34% smaller than the one ob- 
tained by method 1). 

3) The third result, given in rows labeled "combined" 
of Table [VlTl is a x^ fit of a constant to the six R\.j 
values measured in the six proper decay time bins. 
The six R\j values and the result of the fit for each 
of the four subsamples are shown in Fig. [13 The 
statistical uncertainty of this result is 2-3% smaller 
than in method 2), because additional information 
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on proper decay time is included through the six 
eks/eWs ratios. 

Using the MC simulation, we verified that method 3) 
has the best sensitivity; we therefore quote as our final 
result the value obtained by method 3) . To illustrate the 
effect of including the proper decay time information, we 
show also the results of methods 1) and 2). 

From Table IVIII one can see that the central values ob- 
tained by using the three methods, are slightly different. 
In evaluating the significance of the difference between 
the result of method 1) and the result of method 2), we 
have accounted for the ratio of the proper decay time ef- 
ficiencies and for the statistical correlation between the 
samples. The differences are within the expected sta- 
tistical fluctuations: they range between —0.6 and -1-1.4 
standard deviations. By using toy MC simulation, it has 
also been verified that the differences in the central val- 
ues between methods 2) and 3) are within the range of 
expected statistical fluctuations. For the default method 
3), results for Rm in all four subsamples are consistent 
with the null value. The values, shown in Fig. 1151 are 
in good agreement with the expected distribution for 
5 d.o.f., which has a maximum at the value of 3.0. 

The combined result for the electron decay mode is 
obtained by a fit to the values for the e-1 and e-2 
subsamples, obtained by method 3). The flt yields R'j^j = 
(—0.6 ±2.6) X 10""' with a value of 0.1 per 1 degree of 
freedom. The combined result for the muon decay mode 
is obtained in the same way; the fit yields R'^j — 
(5.9 ± 3.7) X 10-'' with a value of 0.4. 

The combined result, taking into account the statistical 
uncertainty only, is obtained by a fit to the four values 
(electron and muon decay mode, SVD-1 and SVD-2); it 
yields a value of 

Rf}"'- = (1.6 ± 2.2) X 10-^ (22) 

where the quoted uncertainty is statistical only. The x^ 
value is 2.5 for three degrees of freedom. The Rm values 
for the four subsamples and the result of the fit are shown 
in the left plot of Fig. [TBI To obtain the final result, the 
partially correlated systematic uncertainties have to be 
studied and taken into account. 



A. Systematic uncertainties 

In the following subsection, different sources of the 
systematic uncertainties are discussed and the system- 
atic uncertainties are given; they are summarized in Ta- 
ble |VlTIl 

1. Finite statistics of the fitting distributions 

One of the main sources of systematic uncertainty is 
the limited statistics of the samples used to obtain the 
signal and background AM distributions used in the 



TABLE VII: The number of fitted signal events in the RS and 
WS samples, the ratio of RS and WS t^y efficiencies, and the 
resulting R\,i value for each proper decay time interval for the 
four subsamples. The results of the fit to the six individual 



R\f values 


are denoted 


as "combined". 




txy 




-^*^WS ^RS/^WS 








^-1 siibsaTTinlp' 




1 fi-9 n 




4:.oziiZii.Zi u.yj-ijziiu.uui 




2.0-2.5 


1 1 97o_|_Qq 




6.1±14.8 








7.2±12.6 


3.1-4.0 


7937±83 


-28.0±25.9 0.310±0.003 


-10.9±10.1 


4.0-5.6 


6394±85 


-21.2±28.6 0.223±0.003 


-7.4±10.0 


5.6-9.0 


4196±89 


15.9±29.8 0.223±0.003 


8.4±15.8 


combined 






-1.7±5.2 


1.6-9.0 


51325±213 


-11.5±65.4 0.413±0.001 


-0.9±5.3 


3II txy 


183496±443 


70.1±141 1 


3.8±7.7 






3-2 subsample: 




1.6-2.0 


32616±150 


-19.1±44.0 0.881±0.003 


-5.2±11.9 


2.0-2.5 


28711±146 


-11.4±41.7 0.603±0.002 


-2.4±8.8 


2.5-3.1 


22513±131 


52.5±41.9 0.415±0.002 


9.7±7.7 


3.1-4.0 


18941±132 


-22.6±41.1 0.285±0.002 


-3.4±6.2 


4.0-5.6 


14796±129 


-18.6±42.3 0.198±0.002 


-2.5±5.7 


5.6-9.0 


9072±128 


25.2±46.5 0.186±0.002 


5.2±9.5 


combined 






-0.1±3.1 


1.6-9.0 


126539±332 


-10.7±102 0.389±0.001 


-0.3±3.1 


all txy 


469947±701 


-369±222 1 


-7.8±4.7 


H-l subsample: 


1.6-2.0 


11314±111 


14.2±34.7 0.921±0.005 


11.6±28.2 


2.0-2.5 


10185±109 


-1.8±33.5 0.637±0.004 


-1.1±21.0 


2.5-3.1 


7893±98 


3.5±30.7 0.440±0.003 


1.9±17.1 


3.1-4.0 


6804±96 


-5.5±31.8 0.303±0.002 


-2.5±14.2 


4.0-5.6 


5350±97 


23.7±33.0 0.214±0.002 


9.5±13.2 


5.6-9.0 


3670±90 


-12.8±35.4 0.217±0.003 


-7.6±20.9 


combined 






2.2±7.1 


1.6-9.0 


45181±245 


-11.2±79.9 0.410±0.001 


-1.0±7.2 


all txy 


163215±485 


-204±180 1 


-12.5±11.0 


/i-2 subsample: 


1.6-2.0 


27612±180 


71.4±54.8 0.876±0.015 


22.7±17.4 


2.0-2.5 


23695±170 


9.3±52.3 0.595±0.010 


2.3±13.1 


2.5-3.1 


18905±154 


82.3±49.8 0.405±0.006 


17.6±10.7 


3.1-4.0 


15488±150 


51.1±50.1 0.273±0.004 


9.0±8.8 


4.0-5.6 


11989±144 


20.4±51.1 0.186±0.007 


3.2±7.9 


5.6-9.0 


7146±138 


-20.3±56.5 0.171±0.016 


-4.9±13.6 


combined 






7.4±4.4 


1.6-9.0 


104556±381 


192±125 0.380±0.002 


7.0±4.5 


txy 


396151±761 


170±284 1 


4.3±7.2 



AM fit to data. To estimate this, we vary the contents 
of all bins of the RS and WS, signal and background 
AAI distributions independently in accordance with each 
bin's statistical uncertainty. We repeat the fit to the RS 
and WS data, calculate the corresponding R\^ in each 
proper decay time interval, and obtain a new Rm value. 
Repeating the procedure 1000 times, the obtained distri- 
bution of Rm values has a Gaussian shape. The sigma of 
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FIG. 13: The AAf distribution of WS events in the six proper decay time intervals for the e-1 (left) and e-2 (right) subsamples. 
The points with error bars are the data, the histogram represents the result of the fit, described in the text, and the small 
contribution on the horizontal axis shows the fitted signal yield. 



the Gaussian, fitted to the distribution, is taken as the 
systematic uncertainty due to the hmited statistics of the 
fitting distributions. The uncertainties are Hsted in Ta- 
ble [Villi hue 1. Larger uncertainties in the muon decay 
mode reflect the fact that, compared to the electron de- 
cay mode, the muon background is larger especially in 
the signal region, and secondly, the embedded slow pion 
sample is smaller due to the smaller Ma value used. Since 
this uncertainty is statistical in nature, it is considered to 
be completely uncorrelated between the four subsamples 
(e-1, e-2, ^-1, 



2. The amount of WS correlated background 

The normalization of the WS-correlated background is 
determined by MC simulation, takings into account the 
central values of branching fractions [23| , of decay modes 
that contribute to this background. 

From MC simulation studies we find that in the elec- 
tron decay mode the largest contributions to the WS 



correlated background come from the following decays: 
K~e^Ve {D^ mesons mainly from Z?*° £)°7 de- 
cays, 34% of the correlated background), K^tt^tt'-' 
(12%), D° ^ K-e+iye7r° (7%), D° ^ K-K+ (7%). In 
total 60% of the correlated background comes from these 
decays. 

In the muon decay mode, the largest contributions 
to the WS correlated background come from 
K-tt+tt'' (19%), ^ K-^i+ly^ (12%), D" ^ 

K-TT+TT^TT^ (7%), D° ^ K~Tr+7T'Tr+ (7%), D° 

K+n-K^ (4%), D° -> K+K-K" (4%). In total 53% 
of the correlated background comes from these decays. 
We calculate the weighted average of the relative uncer- 
tainties of the branching fractions |23|] for the stated de- 
cay modes. For the electron decay mode, the averaged 
relative uncertainty is ±3.6% and for the muon decay 
mode ±5.9%. To take into account the uncertainties of 
the branching fractions used in the MC simulation, we 
repeat the WS fits, changing the amount of the total WS 
correlated background by the average uncertainties on 
the branching fractions. The differences between the re- 
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FIG. 14: Same as Fig. 1131 for the /i-l and /i-2 subsamples. 



suiting Rm values and the default values are taken as the 
systematic uncertainty from this source; they are listed 
in Table Enil line 2. 

This procedure is conservative for two reasons. First, 
by varying the total correlated background instead of 
varying its individual components, the uncertainties on 
the branching fractions are implicitly considered to be 
100% correlated, resulting in the maximum possible sys- 
tematic uncertainty. Second, the modes comprising the 
correlated background contribute significantly also to the 
uncorrelated background. Taking this into account would 
lead to a smaller change in Rm- 

This uncertainty is larger in the muon decay mode, be- 
cause the probabilities to misidentify a pion or kaon as a 
muon are much larger than the corresponding probabili- 
ties for misidentification as an electron. Consequently, in 
the muon decay mode the fraction of the correlated back- 
ground is significantly larger f Table W} . its AM shape 
tends to lower values (Fig. [TT|l and its averaged uncer- 
tainty of the branching fractions is larger. 

The systematic uncertainty from this source is the 



same for the SVD-1 and SVD-2 subsample. Since a signif- 
icant part of the correlated background is due to decays 
common to the electron and the muon decay modes, the 
systematic uncertainties for both decay modes are highly 
correlated. Hence the systematic uncertainty from this 
source will be treated as 100% correlated for all four sub- 
samples (e-1, e-2, /Lt-1, /Lt-2). 



3. The AM shape of the WS uncorrelated background 

We also conservatively account for the uncertainty of 
the AM shape of the WS uncorrelated background. We 
vary Ma within the limits given by the statistical uncer- 
tainties of the ri55 values in Table IVTl The riss statisti- 
cal uncertainty for the data is ±0.18 in the electron decay 
mode and ±0.16 in the muon decay mode. For the em- 
bedded slow pion sample, the riss value changes by 0.1 
for AMa = 5. Hence A/a is varied by ±9 in the electron 
decay mode, and by ±8 in the muon decay mode. With 
the new AAI distributions we repeat the fit to WS data. 
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FIG. 16: Left: the Rm values of the four subsamples with 
the statistical uncertainty only, and the result of the fit to 
these four values (dashed line, x^/d-o.f. = 2.5/3). Right: the 
four R]\f values with the systematic uncertainty included and 
the combined result (dashed line), obtained as described in 
Sec. IV A 71 The dotted lines represent the ±1(t interval. The 
solid line corresponds to no mixing. 



FIG. 15: The resulting R\i values for the four subsamples 
and their average value (dashed line). The dotted lines rep- 
resent the ±1(7 interval. The solid line corresponds to no 
mixing. 

recalculate the Rm values and quote the differences from 
the default values as the systematic uncertainties from 
this source. These uncertainties are listed on line 3 of 
Table mill 

Since A/a is determined for the electron and muon de- 
cay mode separately and the uncertainty on TVa is sta- 
tistical in nature, this systematic uncertainty is consid- 
ered to be completely uncorrelated between both decay 
modes. On the other hand. Ma is determined for SVD-1 
and SVD-2 subsamples together, hence the uncertainty 
is treated as completely correlated between them. 

4- Proper decay time distribution 

To check the reliability of efficiencies ej^g and ratios 
^Rs/^WS' ^^'^ ^'^ estimate the effect of the imperfect fit 
to the proper decay time distribution, the values of e^g 
are compared to an alternative estimate from the fit to 
AM, ej^g^^ = iV^s/iV^°s. This method accounts for the 
influence of the associated signal in the txy distribution. 
In a majority of the txy subintervals, ej^g and Er^^^ typ- 
ically agree within ±2%, the largest discrepancies being 
-9.8% and +4:A%. For the integrated 1.6 < txy < 9.0 
interval, they agree within 0.8%-1.4% for the four sub- 
samples. 

To estimate the effect of the discrepancies, the relative 
difference between e^g and e^^^^ is assigned as the rel- 
ative uncertainty on Ers/^ws- Hence we reduce the six 
effciency ratios simultaneously by this uncertainty and 
repeat the Rm calculation; we then increase the ratios by 
this uncertainty, and again recalculate Rm- The differ- 
ence between the resulting Rm value and the default fit 
is quoted as the systematic uncertainty from this source. 
It is very small and can be found on line 4 of Table IVTIIl 



5. The amount of the associated signal 

The systematic error due to the uncertainty in the asso- 
ciated signal fraction is estimated by varying the fraction 
and repeating the fitting procedure. Taking into account 
the uncertainties on the measured branching fractions 
[1^ of the associated signal decay channels, we conser- 
vatively vary the amount of associated signal by ±40%. 
We recalculate the Rm values and compare them to the 
default Rm value; we quote the differences as the sys- 
tematic uncertainty from this source. From Table IVIIII 
(line 5) one can see that it is almost negligible. 



6. The amount of the RS correlated background 

From MC simulation studies we find that in the elec- 
tron decay mode the largest contributions to the RS 
correlated background come from the following decays: 
/iT^Tr+Tr" (33% of the RS correlated background), 
D° ^ K-TT+Tr°Tr° (14%), D° ^ R-e+i^e (13%). In total 
60% of the RS correlated background comes from these 
three decays. 

In the muon decay mode, the largest contributions to 
the RS correlated background come from D'^ K^tt^tt^ 

(43%), D° ^ K-TT+TT°TI° (17%), ^ K-TI+TT-TT+ 

(12%). In total 72% of the correlated background comes 
from these decays. We calculate the weighted average of 
the relative uncertainties of the branching fractions [2^ 
for the stated decay modes. For the electron decay mode 
the averaged relative uncertainty is ±4.3% and for the 
muon decay mode ±4.4%. 

We repeat the RS fits, changing the amount of the 
total correlated background by the average uncertainties 
on the branching fractions. The differences between the 
obtained values of Rm and the default values are taken 
as the systematic uncertainty from this source. They can 
be found in Table IVini (line 6) and are negligible. 
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TABLE VIII: A summary of the systematic uncertainties on 
Rm, the total systematic uncertainty in each subsample and 
the combined (summed in quadrature) statistical and system- 
atic uncertainty. Values are given in units of 10^^. 



source 


e-1 


e-2 


/i-1 




1 fitting histo. statistics 


±1.54 


±0.91 


±2.64 


±1.81 


2 WS correlated bkg. 


±0.37 


+0.39 
-0.38 


+2.98 
-2.89 


+3.05 
-2.97 


3 WS uncorrelated bkg. 


+1.30 
-1.88 


+ 1.70 
-1.85 


+2.58 
-2.82 


+1.57 
-3.20 


4 imperfect t^cy 


±0.05 


±0.02 


±0.01 


+0.25 
-0.33 


5 associated signal 


+0.01 
-0.00 


±0.00 


+0.09 
-0.10 


±0.02 


6 RS correlated bkg. 


±0.00 


±0.00 


±0.01 


±0.04 


total systematic 


+2.05 
-2.46 


+ 1.97 
-2.10 


+4.75 
-4.83 


+3.89 
-4.74 


statistical + systematic 


+ 5.58 
-5.74 


+3.66 
-3.73 


+8.53 
-8.57 


+5.86 
-6.45 



7. Total systematic uncertainty and the final result 

The final result of the measurement is obtained by av- 
eraging the results for the four subsamples, e-1, e-2, /i— 1 
and fj, — 2. As explained at the beginning of Sec. |Vl the 
results obtained by method 3) are used (quoted in Table 
NTH as "combined". 

The contributions to the systematic uncertainty are 
divided into three categories: 

(a) The systematic uncertainty that is completely cor- 
related between all four subsamples. This is the 
error due to the uncertainty of the WS correlated 
background fraction. Sec. IV A 21 

(b) The systematic uncertainty that is completely cor- 
related between the SVD-1 and SVD-2 subsamples 
and is uncorrelated between the electron and the 
muon decay mode. Such a contribution comes from 
the uncertainty of the AM shape of the uncorre- 
lated WS background. Sec. IV A 31 

(c) The systematic uncertainties that are uncorrelated 
between the four subsamples, or are very small. 
The main contribution comes from the uncertainty 
due to the finite statistics of the fitting distribu- 
tions. Sec. IV A II The uncertainties from all the 
remaining sources are also added. 

To obtain the final result and its uncertainty, taking into 
account the systematic uncertainties, we adopt the fol- 
lowing procedure: 

(1) For each of the four subsamples, we add to the sta- 
tistical uncertainty in quadrature all the uncertain- 
ties from category (c). 

(2) We perform the fit to the SVD-1 and SVD-2 
Rm values in the electron and muon decay mode 
to obtain the averaged value for the electron and 



muon decay mode, (-0.56 ± 2.76) x lO""* and 
(5.89 ±4.02) X lO"-*, respectively The quoted un- 
certainties include the statistical uncertainty and 
the uncertainties (c). 

(3) To add the uncertainty (b) for the electron decay 
mode, we first simultaneously increase and then si- 
multaneously decrease the results for the e-1 and 
e-2 subsamples by the uncertainty (b) and repeat 
step (2). The difference from the default result 
of step (2) is added in quadrature to the uncer- 
tainty obtained in step (2). The result for the elec- 
tron decay mode, including the statistical uncer- 
tainty and systematic uncertainties (b) and (c), is 
(— 0.56t3;33) X 10"''. We perform the same proce- 
dure also for the muon decay mode; the result is 

(5.891^^0?) X 10-^. 

(4) We perform a fit to the results for the electron 
and muon decay mode, obtained in step (3); the 
result is (1.27 ± 2.70) x 10"^. The obtained mean 
value is the final result, but the uncertainty needs 
to be increased by the uncertainty (a). 

(5) To account for the uncertainty (a) we first simulta- 
neously increase and then simultaneously decrease 
the initial four Rm values by the uncertainties (a) 
and repeat the steps (l)-(4). The difference from 
the default result, x 10~^, is added in quadra- 
ture to the uncertainty previously obtained from 
step (4) to obtain the final uncertainty of the re- 
sult. 

The total uncertainty of the final result is ±2.93 x 10~^. 
We calculate the contribution of the systematic uncer- 
tainty as the difference between the total uncertainty and 
the statistical uncertainty (Eq. 2.93^ - 2.16^ = 

1.98^. The final result is then 

= (l-3±2.2±2.0) X 10-^ (23) 

where the first uncertainty is statistical and the second 
systematic. As this value is close to the boundary of the 
physical region {Rm > 0) we use the Feldman-Cousins 
approach [26|] to calculate upper limits: 

Rm < 6.1 X 10"^ at the 90% confidence level, (24) 
Rm < 7.0 X 10"^ at the 95% confidence level. (25) 

With systematic uncertainties included, the final results 
for the electron and muon decay modes are: 

= (-0.6 ± 2.711;?) X 10-^ (26) 
i?^, = (5.9 ±3.7l3 J) X 10-4. (27) 

The Rm values for the four subsamples, including the sys- 
tematic uncertainty, and the combined result are shown 
in the right plot in Fig. [161 

The increase in the sensitivity of the current result, 
compared to the one published in 0], is caused partially 
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by the larger statistical power of the sample, but also 
by the improvements in the measurement method. The 
statistical uncertainty of the present result in the elec- 
tron sample (Eq. ([26]) ) is about 22% smaller than one 
would expect by appropriately rescaling the uncertainty 
of the result [7] by the increase of the data set used. The 
improvement is mainly due to improved selection crite- 
ria, improved neutrino reconstruction and improvements 
in using the D'^ proper decay time measurement. The 
systematic error of the result in the electron sample is, 
however, larger than the one published in [rj as it is es- 
timated more conservatively. 



VI. SUMMARY 

Using a data sample with an integrated luminosity 
of 492.2 fb~^, collected by the Belle detector, we have 
searched for D^-D^ mixing using semileptonic decays 
of the neutral charmed meson, DO K'^*)+e-Ve and 
Z?" K'^*^^ fi^Pfj^. We select mesons produced via 
the decay — > tt^Z?", and tag the flavor of the D 
meson at production by the charge of the accompanying 
slow pion. The measured mixing rate Rm is consistent 
with no mixing in both electron and muon decay modes. 
The combined result accounts for the partially correlated 
systematic error and yields Rm = (1.3±2.2±2.0) x 10"''. 
Since it is consistent with zero we set upper limits on the 
mixing rate of Rm < 6.1 x 10~^ at the 90% confidence 
level. 

This result supersedes that published in Ref. 7] and 
represents the most stringent experimental limit on R]\j 
obtained to date from semileptonic decays. Its ac- 
curacy is significantly better than that of the world av- 
erage of previous measurements in semileptonic decays, 
Rm = (1.7 ±3.9) X 10"'' Q. Although the sensitivity 
is not sufficient to observe a positive mixing signal, it is 
worth noting that in semileptonic decays no model un- 
certainties can influence the result. The reported value 
of Rm is in agreement with the world average values of 
X = (0.87t°j^)% and y = (0.66t°JJ)% [| and it wiU 
help in further constraining the D'^ mixing parameters 
in combination with the results of the measurements in 
other decay channels. 

Appendix 



In the proper decay time fit the following functions are 
used: 

The Lorentz function, centered at 0: 



L{t;b)^- ^ 



TT 1 + (te)2 

The asymmetric Lorentz function, centered at 0: 



(28) 
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